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W). Abstract 

We re- formulate the 3 + 1 GRMHD equations for the Schwarzschild 
j> ■ black hole in a Veselago medium. Linear perturbation in rotating 

(non-magnetized and magnetized) plasma is introduced and their Fourier 
_J | analysis is considered. We discuss wave properties with the help of 

£S) ■ wave vector, refractive index and change in refractive index in the 

r^| . form of graphs. It is concluded that some waves move away from 

the event horizon in this unusual medium. We conclude that for the 
rotating non-magnetized plasma, our results confirm the presence of 
Veselago medium while the rotating magnetized plasma does not pro- 
vide any evidence for this medium. 
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1 Introduction 

Schwarzschild black hole is the simplest one which has neither angular mo- 
mentum nor charge. Plasma is the most common form of matter in the 
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universe. It can be accelerated and directed by electromagnetic fields which 
make it to be controlled and employed. The matter in inner-stellar space is 
composed of plasma. The general relativistic magnetohydrodynamics (GRMHD) 
is the standard theory which is helpful to explore the dynamics of falling 
magnetized plasma towards the event horizon. 

An accretion disk is formed near the Schwarzschild black hole as plasma 
flows towards it. Perturbations in black hole regime have always been fas- 
cinated the researchers. Regge and Wheeler [1] discussed the stability of 
Schwarzschild singularity. Price [2] examined the dynamics of non-spherical 
perturbations during the collapse of stars with a scalar-field analog. Waves 
in an electron-positron plasma were investigated by Sakai and Kawata [3] in 
a frame of two-fluid equation for the Schwarzschild black hole. Southwood 
[I] used a small perturbation approach to derive the general stability criteria 
for plasma with isotropic pressure. Gleiser et al. [5] studied even parity 
perturbation of the Schwarzschild black hole up to second order. 

Arnowitt, Deser and Misner (ADM) [6] developed 3+1 formalism in which 
spacetime is foliated into layers of three-dimensional spacelike hypersurface, 
threaded by a timelike normal. Many relativists [TJ- [5] applied this technique 
to illustrate different attributes of general relativity. Thorne and Macdonald 
[9]- [10] evolved the electromagnetic fields of the black hole theory. Holcomb 
and Tajima [UJ, Holcomb [12] and Dettmann et al. [13] analyzed the pro- 



duction of waves in the Friedmann universe. Buzzi et al. [H] examined 
the wave properties of two fluid plasma in the surroundings of Schwarzschild 
event horizon. Zhang [15] explicated the theory of stationary symmetric 
GRMHD in black hole regime. The same author [16] also studied the behav- 
ior of perturbation of cold plasma in the magnetosphere of Kerr black hole. 
Sharif and his collaborators [I7]-[2U] determined the plasma wave properties 
by using dispersion relations. They analyzed it for cold, isothermal and hot 
plasmas. 

Veselago medium or Double negative medium (DNG) is the most famous 
class of metamaterials. This medium has simultaneously negative electric 
permittivity and magnetic permeability. Negative phase velocity (NPV) and 
negative refractive index are main features of this medium. Ziolkowski [2T] 
analyzed electromagnetic wave properties in DNG. Valanju et al. [22] estab- 
lished the well-known idea that wave refraction in DNG medium is always 
positive but not homogenous. Ross et al. [23] found that NPV propaga- 
tion develops at lower values of cosmological constant A in the vicinity of 
non-rotating black hole. Mackey and Lakhtakia [24] verified that presence 



of charge raises the tendency of a rotating black hole to support the NPV 
propagation in its ergosphere. In a paper given in a book [25], J li general- 
ized the concept of negative medium to acoustic waves. In a recent paper, 
Veselago [26] has explained how an electromagnetic wave transfers energy, 
linear momentum and mass in a negative refraction-medium. 

In the present paper, we discuss isothermal plasma wave properties of 
the Schwarzschild black hole in a Veselago medium. The format of the paper 
is as follows: General line element and its formation in the Schwarzschild 
planar analogue is given in Section 2. In Section 3, perturbed and Fourier 
analyzed forms of 3+1 GRMHD equations for isothermal plasma are pre- 
sented. Sections 4 and 5 contain the restricted 3 + 1 GRMHD equations for 
rotating (non-magnetized and magnetized) plasmas and discussion related to 
wave properties. Summary of the results is given in Section 6. 

2 ADM 3+1 Formalism and Schwarzschild 
Planar Analogue 

The general line element in ADM 3+1 formalism is given by [16] 

ds 2 = -a 2 dt 2 + j ij {dx i + (3 i dt)(dx j + ftdt), (2.1) 

where the lapse function a is the ratio of FIDO (fiducial observer) proper time 
to universal time, i.e., a = %, (3 l is the shift vector and 7^ (i,j = 1, 2, 3) are 
the components of three dimensional absolute space. FIDO is the natural 
observer associated with the above line element. In Schwarzschild planar 
analogue, shift vector vanishes due to zero angular momentum and ja = 1 
represents Euclidean 3 geometry. Thus it reduces to the form [T7] 

ds 2 = -a 2 (z)dt 2 + dx 2 + dy 2 + dz 2 , (2.2) 

where the directions z, x and y are analogous to the Schwarzschild coordi- 
nates r, (f> and 9 respectively. 

3 3+1 GRMHD Equations for Isothermal Plasma 
in Veselago Medium 

The 3 + 1 GRMHD equations in a Veselago medium for the plasma present 
in general line element and the Schwarzschild planar analogue (Eqs. (j2.ip and 



( 12. 2p ) are given in Appendix A. In the vicinity of the Schwarzschild black 
hole, equation of state for isothermal plasma is 



(i= , (3.1) 

Po 

here po, p, P and \i are the rest mass density, moving mass density, pressure 
and specific enthalpy respectively. The specific enthalpy p is constant here. 
This equation indicates that there is no exchange of energy between the 
plasma and the magnetic field of the fluid. For isothermal plasma existing 
in Schwarzschild magnetosphere, the 3+1 GRMHD equations (A10)-(A14) 
take the form 

_ = -V x (qV x B), (3.2) 

V.B = 0, (3.3) 
^p) + («V.V)(p + p) + (p + p) 7 2 V.^ + (p + ph 2 

V.(aV.V)V + (p + p)V. (aV) = 0, (3.4) 

{ (<" + ^ + S) s -> + { » + >* v < v > ~ h 3 *) 61 

+V.V) V - H^Sy - ^B r B,^\ ViV k + l> + p)j 2 ai + p,, 

= ^(V x B),V.(V x B) - JLjCaB); + ^(aflO ^ 

--L[B x {V x (V x (aV x B))}] ( , (3.5) 

Aira 

(-^ + V - V )^ + ^ ~ ~ tr + 2 0° + P)7 2 (V.a) + (p + p) 
a at a eft 

1 r)R 1 r)B 

T'(V.V) - ^ (V x B).(V x _) - — (V x B).(B x _) 

+- — (VxaB) = 0. (3.6) 

Alia 

We assume that plasma flows in xz-plane so that the velocity V and 
magnetic field B experienced by FIDO are given as 

V= V(z)e x + u(z)e z , 

B = B[X(z)e x + e z ], (3.7) 



where B is an arbitrary constant. The relation among A, u and V is given 

by mi 

V F 
V = — + Xu, (3.8) 

a 

where V F is a constant of integration. 7 = , 1 is the Lorentz factor 
which becomes 

The gravity of black hole perturbs the plasma flow. Inserting the linear 
perturbation to density p, pressure p, velocity V and magnetic field B, we 
have 

p = p° + 5p = p° + pp, p = p° + 5p = p° + pp, 

V = V° + <5V = V° + v, B = B + <5B = B° + J Bb, (3.10) 

where unperturbed and linearly perturbed quantities are denoted by p°, p, V°, 
B° and 5p, 6p, <5V, 5B respectively. The quantities p, p,v x ,v z ,b x and b z are 
dimensionless which we shall introduce for the perturbed quantities 

p = p(t,z), p = p(t J z) J 

v = 5V = v x (t, z)e x + v z (t, z)e z , 

b = — = b x (t,z)e x + b z (t 1 z)e z , (3.11) 

The perfect GRMHD equations (Eqs. fl3.2l) - f)3.6p ). after the insertion of linear 
perturbation from Eq. (13.111) . become 

-ip = -Vx(ftvxB)-Vx(«Vx5B), (3.12) 

V.(5B) = 0, (3.13) 

d -^P±^L + {a\.V){5p + 5p) + {p + p) 1 2 \.^ [ -a{p + p) 

(v. V In u) + a(p + p) ( V. v) + (5p + Sp) (V.aV) + (5p + 5p)^ 2 
V.(aV.V)V + 2(p + p) 7 2 (V.v)(aV.V) In 7 + (p + p)-f 2 
(aV.VV).v + (p + p)7 2 V.(aV.V)v = 0, (3.14) 



{((p + ,)7 2 + g^ + (p + p)7v^.-^,}i^ 

+i-[B x {V x ^^}], + (p + ptfvuV* + (p + P ) 
x^ Vj , k V j V k - -L{(a6Bi)j - (a5B,)^}B ] + (5p) t 

47TQ; 

= - 7 2 {(5p + 5p) + 2(p + p) 7 2 (V.v)H + ^-{(aB^ 

-(aBj) ti }8BS -(p + p)j 4 ( Vi V j + viVi)V kj V k - 7 2 {(<5p + Sp) 
V 1 + 2(p + ph 2 (V.v)V j + (p + p)v j }V itj - i 4 Vi{(5p + Sp)V j 
+4(p + p) 7 2 (V.v)^' + (p + py'}V,-*V*, (3.15) 



2 ld(5p + 6p) , _ , , 2 10(<Jp) 



a, at a at 



+2(p + p) 7 4 (V.V)(V.v) + 2(p + p) 7 2 (v.a) + 4(p + p) 7 4 (V.v) 

(V.a) + 2(5 p + <5p) 7 2 (V.a) + (p + p) 7 2 (V.v) + 2(p + p) 7 4 

1 r)B r)B 

(V.v)(V.V) + (Sp + 5p)7 2 (V.V) = _ [v.(B.— )V + V.(B. — )v 

f)fi f)B /MB 

+V.(B.5B)V + V.(5B—)\ - v.(B.V)^- - V.(B.V)^- 

_V,B.v)^-V,,B.V)f 1 -J, [ V,B.B)|-V,B.> 1 

— V x (aSB). (3.16) 

Ana 

The component form of these equations, using Eg. ( 13. 11]) . can be composed 
as follows 

— xr ~ ub x,z = (ub x ~ Vb z -v x + Xv z )V In a 
a at 

~(v x , z - Xv ZyZ - X'v z + V% + Vb z , z - u%), (3.17) 

1 db z , 

-of = 0, 3.18 

a at 

b z , z = 0, (3.19) 



(JO (JT) (JV (JV 

P ~dt + P ~dt +iyp + P ^ 2(yV ~dt~ + U ~dt~" > + aupp > z + au W,z 

+a(p + p){rfuVv X)Z + (1 + 7 2 m 2 )w 2 , 2 } (p-p)(awyp) tZ 
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+a{p + p)i 2 u{{l + 2 1 2 V 2 )V + 2 1 2 uVu'}v x - a{p + p) 

v' 
x{(l - 2 7 V)(1 + 7 V)-i - 2 7 4 uVy> 2 = 0, (3.20) 

u 

( (p+ph2(1+7 ^ ) + ^}I^ + ( (p+ph w-^ 

X a1^ + {(P + Ph 2 (l + 7V 2 ) + f4 w*,* + {(P + P)7 4 ^ 

\ 7^2 "j 7^2 7^2 

m> z , 2 - — (1 + u 2 )b XtZ - - — {a'(l + u 2 ) + aW} 6^ 



An J ' 47r ' 47ra 

+7 2 M(pp + pp) { (1 + 7 V 2 )\/' + 7 2 wlV} + [(p + p) 7 4 m{(1 

B 2 W 

+4 7 2 \/ 2 )W + AVV'(1 + ! 2 V 2 )} + }v x + [(p + p) 7 2 {(l 

Ana 

+2 7 V)(1 + 2 1 2 V 2 )V - 7W + 2 7 2 (1 + 2 7 V)wW} 

£ 2 ^ xn , , 

f, x 2,, 2 2^ A 2 £ 2 1 19U* f, x 4 t, ^S 2 

x a if + 1 (p + p ' 72 ' 1 + 72 " 2 ' + ^nr 1 "" w + " p + p ' 74 " 

xV- -;—}««*,* + -;— (1 + «)&*,* + - — {a X - (aX)' + wA 

47T 47T Alia 

x (ua' + u'a)}^ + (pp + Pp)l 2 {a z + «w'(l + 7V) + 7V W'} 
+ [(p + p)7 4 {mV(1 + 4 7 2 \/ 2 ) + 2\/(a, + W(l + 2 7 V))} - XB 2 

x T^K + [(P + P)7 2 {«'(1 + 7 V) (1 + 4 7 V) + 2« 7 2 (a 2 + (1 
47ra 

+2 7 V)VV')} + ^(«A)> 2 + (p'p + pp') = 0, (3.22) 



1 2 dp 1 2 dp 2 . . ., 2 . _ _ . 

~7 P a7 + ~7 P aT + 7 IP + P K + M 7 (PP,z + PP,z + PP 
a at a at 

1 OX) 

+p'p) P^rr + 2l 2 u(pp + pp)a z + j 2 u'(pp + pp) + 2(p 

a at 

+p)j i {uV + 2uVa z + u'V)v x + 2(p + p)7 2 (2 7 2 W + a z7 4 
+2 7 Va> 2 + 2(p + p)i A uVv x>z + (p + p) 7 2 (l + 2 7 V) 

xv v - - — [(V 2 + u 2 )\b x + (V 2 + u 2 )b z - \V(\V 

+u)^ - u(XV + u) 8 ^} - ^[(V - \u)v x>t + \{u\ 

-V)v Zjt ] + ^-b x , z = 0. (3.23) 

For Fourier analysis, we take the harmonic spacetime dependence of pertur- 
bation 



p(t, z) 


= cie- t( ^- fcz) , 


p{t,z) = c 2 e~ i{uJt - kz \ 


v z (t,z) 


= c 3 e~ L{uJt - kz) , 


v x (t,z) = c,e-^ t - kz ^ 


b z (t,z) 


= c 5 e- L{uJt - kz) , 


b x (t,z)=c 6 e-^ t ^ 



(3.24) 

Here k is the z-component of the wave vector (0, 0, k) and u is the angular 
frequency. Using the wave vector, we obtain refractive index which helps to 
examine the behavior of plasma waves near the event horizon. Wave vector 
and dispersion relation can be defined as 

• Wave Vector: A vector whose direction indicates the direction of 
phase propagation of a wave is called wave vector. Its magnitude is the 
wave number. 

• Dispersion Relation: Dispersion relation gives the angular frequency 
u as a function of wave vector k in the form ^ = Xu, where A is the 
wavelength [27]. Dispersion is said to be anomalous if the refractive 
index is less than one and its change with respect to angular frequency 
is negative, otherwise normal. 



The Fourier analyzed form of Eqs. fl3.17p - fl3.23p . by using Eq.f l3.24p . are 

c±{a' + ika) — c 3 {(a\)' + tka\ } — c 5 (aV)' + c 6 {(au)' + uo 
+ikua} = 0, (3.25) 

c*(— ) = 0, (3.26) 

a 

c 5 ik = 0, (3.27) 

ci{(—ioj + ikau)p — p^ 2 au{VV' + uu') — a' up — au'p — aup'} 
+C2{(— uv + ikau)p + a' up + au'p + aup' + pj 2 au(VV' + uu')} 
+c 3 (p + p) [-LL0-y 2 u + ika{l + 7V) - a{(l - 2j 2 u 2 )(l + -y 2 u 2 ) 

u' 

x 27VW}] + c 4 (p + p)h 2 V{-Luo + ikau) + «7 2 m{(1 

u 

+2 1 2 V 2 )V + 2 1 2 uVu'}] = 0, (3.28) 

c lP1 2 u{(l + 7 V 2 )F' + ~i 2 uVu'} + c 2 p 7 2 M{(l + 7 V 2 )y' 

\ d2 .. . 

+7 2 mVm'} + c 3 [-{(p + p) 7 4 m1/ — }— + Up + p)7 4 mV" 

47T a 

-^-}^w + (p + p)7 2 {(l + 27V) (1 + 2 7 V 2 ) - 7 V 2 }V 
An 

B 2 u 
+2 7 4 (p + p)uFu'(l + 2 7 V) - - — (aX)'} + c 4 [-{{p + j9) 7 2 

Ana 

(1 + 7V 2 ) + ^}- + {(p + Ph 2 (l + 7V 2 ) + ^}iku 

An a 47T 

£? 2 7/r/ 

+ (p + p) 7 4 w{(l + 4 7 V 2 )W + AVV'(1 + 7V 2 )} + ] 

Ana 

B 2 a' 

-c 6 — {(1 + u 2 )ik + (1 + u 2 )— + W} = 0, (3.29) 

An a 

c\pl 2 {a z + uu'(l + 7 V) + 7VW} + c 2 [p7 2 {a 2 + W 

(1 + 7 V) + 7VW} +p' + ikp] + c 3 [-{(p + p) 7 2 (l + 7V) 

\2d2 A 2 5 2 

+ -—}— + {(p + P)l 2 (l + l 2 u 2 ) + ——}iku + {(p + P h 2 
An a An 

x{u'(l + 7V)(1 + 4 7 V) + 2w 7 2 {a, + (1 + 2j 2 u 2 )}W} 
+ ^(«A)'] + c 4 [-{(p + ph%V - ^}^ + {{p + ph%V 



XB 2 

" -}iku + {{p + ph 4 {u V(l + 4 7 V 2 ) + 2V{(1 + 2 7 V)m/ 



47T 

\ d2„/„ C?2 

+a z — - — }] + c 6 [— — {—(a A)' + a'A - uX(ua' + w'a)} 

411a 411o; 

\d2 

+— (l + u 2 )^] = 0, (3.30) 

c i{(— ^7 2 + ^«7 2 + 2u7 2 a 2 + 7 V)p + Vt 2 } + c 2 {( — (1 - 7 2 ) 
a a 

+tku^f 2 + 2^ 2 ua z + 7 2 M ; )p + wy 2 p'} + c 3 7 2 {(p' + p') + 2 

x (2 7 W + a 2 + 2 7 Va 2 )(p + p) + (1 + 2 7 V)(p + p)ifc + 

Ana 

x (Am - V)lu] + c 4 [2(p + p) 7 2 {(uK' + 2uVa z + uV) + uVtk} 

B 2 —B 2 
+- — (V - u\)uj] + ce[- {(^ 2 + w 2 )A + XV (XV + u)iw} 

+^Lk] = 0. (3.31) 

These equations will be used to find dispersion relations. 

4 Plasma Flow with Rotating Non- Magnetized 
Background 

In the present section we shall study the rotating non-magnetized plasma, 
i.e., B = 0. Equations (13.21) and (13.31) which are evolution equations of the 
magnetic field vanish. In the Fourier analyzed perturbed GRMHD equations 
(Eqs. (13.28jl - (13.3ip ). we put B = = A and C5 = = cq and obtain 

ci(— p) + c 2 (—p) + c 3 (p + p)[— 7 2 « + (1 + J 2 u 2 )tk 
a a a 

-(1 - 2 7 V)(1 + 7 V)- + 2 7 4 «W] + c 4 (p + ph 2 [(— 

u a 

+iku)V + u(l + 2 7 V 2 )V + 2 7 2 uVm'] = 0, (4.1) 

ciP7 2 m{(1 + -i 2 V 2 )V + 7 2 VW} + c 2 P7 2 m{(1 + ^V 2 )V 

+>y 2 Vuu'} + c 3 (p + p)7 2 {(— — + Lku)-y 2 Vu + (1 + 2 7 V) 

a 

x(l + 2 7 V 2 )V" - 7 W + 2 7 2 (1 + 2 7 V)uW} + c 4 {p 
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+ph 2 [(l + 1 2 V 2 )(— — + iku) + 7 2 m{(1 + 4 7 VW + 41/ 
a 

V'(1 + 1 2 V 2 )}} = 0, (4.2) 

CiP7 2 K + (1 + tV)W + l 2 u 2 VV'} + c 2 [p-i 2 {a z + (1 + 7 V) 

xW + tVW} + (p' + tfcp)] + c 3 (p + p)7 2 [(l + 7V)(— — 

+i£;w) + u'(l + 7 V) (1 + 4 7 V) + 2«7 2 {a z + (1 + 2 7 V) 

x\/\/'}] + c 4 (p + p)^\uV(— + Lku) + u 2 V'(l + 4 7 V) 

a 

+2V{a z + (1 + 2t 2 u 2 )mw / }] = 0, (4.3) 

c i{(— — 7 2 + tAwy 2 + 2u7 2 a z + 7 2 «')P + Vt 2 } + c 2{( — (1 - 7 2 ) 

+iku>y 2 + 27 2 wa 2 + 7 2 w')p + «7 2 p'} + c 3 7 2 {(p' + p') + 2(2j*uu' 
+a z + 2 7 Va 2 )(p + p) + (1 + 2 7 V) (p + p)ik} + c 4 {p + p)2 7 4 
{{uV + 2uVa z + u'V)+uVik} = 0. (4.4) 

4.1 Numerical Solutions 

In order to determine numerical solutions, we use the the following assump- 
tions 

• Specific enthalpy: p = 1. 

• Time lapse: a = tanh(10^)/10, 

• Stationary fluid: aj = 1 with velocity components V — u yields the 
following expression: aj = 1 => 7 = 1/Vl — u 2 — V 2 = l/a. 

• Velocity components: u — V, x and z-components of velocity lead to 
U = V=- ^~" 2 



With the assumption of stiff fluid, i.e., p = p, when we replace these val- 
ues, the mass conservation law in three dimensions yields p = p = —l/2u. 
The GRMHD equations (Eqs. (j3.2j) - (13.6p ) are satisfied by the above assump- 
tions for the region 1.5 < z < 10, < 00 < 10. For plasma flow in ro- 
tating non-magnetized background, we acquire the constant values of u — 
V = —0.703562 which make the flow constants I = 0.703562 and e = 1. 
The determinant of the coefficients of constants of Eqs . (14 . 1 [) - (14 . 41) are solved 
which lead to a complex dispersion relation [28]. By comparison of real and 

11 



imaginary parts, two dispersion relations are found. The real part of the 
determinant gives an equation quartic in k 

A^z)^ + A 2 (z, u)k 3 + A 3 (z, u)k 2 + A A (z, u)k + A 5 (z, w) = (4.5) 

which gives two real values of k and two complex values conjugate. The 
equation obtained from the complex part is cubic in k 



B x {z)k? + B 2 (z, u)k 2 + B 3 (z, u)k + B 4 (z, u) = 



(4.6) 



which yields one real value of k and the remaining two are complex conjugate 
of each other. We can compute refractive index and its change with respect 
to angular frequency by using the real values obtained from Eq. fl4.5|) and 
( 14.61) . The results reached are shown in Figures 1, 2 and 3. The following 
table shows the results obtained from these Figures. 

Table I. Direction and refractive index of waves. 





Direction of Waves 


Refractive Index (n) 


1 


Move towards the event horizon 


n < 1 in the region 

< z< 10,8.5 <lu < 10 and 

0< z < 10,9.5 <lu < 10 

with the decrease in z 


2 


Move towards the event horizon 


n < 1 in the region 

2 < z < 10, 2 < u < 9 and 

2< z < 10,9.9 <u < 10 

with the decrease in z 


3 


Move away from the event horizon 


n < 1 in the region 
0< z < 10,2 <u < 10 
with the decrease in z 



In Figure 1, the change in refractive index with respect to angular fre- 
quency represents normal as well as anomalous dispersion of waves at random 
points. In Figure 2, dispersion is anomalous in the regions 2 < z < 10, 2.2 < 
uj < 2.25 and 0<z<9,7.5<u;<8 while Figure 3 also shows normal and 
anomalous dispersion randomly. The group and phase velocities are found 
to be antiparallel in all figures. 
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Surface Plot of Wave Vector 




Surface Plot of Phase Velocity 



Surface Plot of Group Velocity 



0.0703562 

0.0703562' 

0.0703562' 

0.0703562 

0.070356" 




-0.0703562 
-0.0703562' 
-0.0703562' 




Surface Plot of Refractive Inde. 




Figure 1: Waves move towards the event horizon. The dispersion is normal 
as well as anomalous at random points. 
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Surface Plot of Wave Vector 




Surface Plot of Phase Velocity 



0.0703562 

0.0703562' 

0.0703562 1 

0.0703562' 

0.0703562' 




Surface Plot of Group Velocity 



-0.0703562' 

10 -0.0703562' 
-0.0703562' 
-0.0703563 1 




Surface Plot of Refractive Index 



Surface Plot of D[n,(j] 




10 -0.00001 
-0.00002' 
-0.00003 




Figure 2: Waves are directed towards the event horizon. Region has anoma- 
lous dispersion. 
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Surface Plot of Wave Vector 




Surface Plot of Group Velocity 



Surface Plot of Phase Velocity 




-0.109255 1 
-0.109255' 
-0.109255' 
-0.109255' 




Surface Plot of Refractive Index 




Figure 3: Waves are directed away from the event horizon. The region has 
normal and anomalous dispersion of waves. 
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5 Plasma Flow with Rotating Magnetized Back- 
ground 

This is the general case in which we assume that plasma is magnetized and 
rotating. We suppose xz-plane for the velocity and magnetic field of fluid. 
Equations (I3.25I) - (I3.31I) are the corresponding Fourier analyzed perturbed 
GRMHD equations. 

5.1 Numerical Solutions 

For magnetized background, we assume V F = with u = V so that Eq. (13.81) 
gives A = 1. Also, we consider B = y^y- which leads to ^ = 2. Further, 
the values of lapse function, velocity, pressure, density and specific enthalpy 
remain the same as given in section 4.1. The above assumed values satisfy 
the perfect GRMHD Eqs. <ET2)l - (EP)]) for the range 1.5 < z < 10, < u < 10. 
The values of flow constants are / = e = 1, I — —7.32001 and h = —1.76321. 
Also, we have u — V — —0.703562. We obtain the determinant of the 
coefficients of constants by substituting these values in Eqs. ()3.25p and (13. 28ft - 
(I3.3ip which gives two dispersion relations. Also, C5 = from Eqs. (l3.26p - 
( I3.27p . The dispersion relation obtained from the real part has the following 
form 

Ai(»£; 4 + A 2 (z, u)k 3 + A 3 (z, u)k 2 + A A (z, u)k + A 5 (z, w) = (5.1) 

while the imaginary part gives the following dispersion relation 

B x (z)k 5 + B 2 {z, u)k A + B 3 {z, u)k 3 + B 4 {z, co)k 2 + B 5 (z, u)k 
+B 6 (z,uj) = 0. (5.2) 

We have used the software Mathematica to solve the above equations so that 
the roots can be displayed in terms of graphs. However, we could not find 
any graph of the roots. It seems that either all the roots are imaginary or 
there does not exist any wave in this region. 

6 Summary 

This paper investigates the isothermal plasma wave properties of the Schwarzschild 
magnetosphere in a Veselago medium. For this purpose, we have re-formulated 
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the Maxwell and 3 + 1 GRMHD equations by considering both permittivity 
and permeability less than zero. The component and Fourier analysi of these 
equations are derived for rotating plasma in non-magnetized and magnetized 
backgrounds. 

For the rotating non-magnetized plasma, our assumed values satisfy the 
3 + 1 GRMHD equations in the region 1.5 < z < 10. In Figures 1 and 2, 
waves move towards the event horizon. The dispersion is normal as well as 
anomalous at random points in the first figure while it is anomalous in the 
second figure. The third figure indicates that waves move towards the event 
horizon while it has normal and anomalous dispersion of waves randomly. All 
the figures show that phase and group velocities are antiparallel. Refractive 
index is less than one and it increases in a small region near the event horizon. 
From the previous literature [19] of the usual medium for isothermal plasma, 
we know that all waves move towards the event horizon. In the Veselago 
medium, Figure 3 indicates that waves move away from the event horizon. 
This is entirely a different result in this unusual medium. For the usual 
medium, the value of refractive index is always greater than one while it is 
less than one in all figures. This confirms the presence of Veselago medium. 
The rotating magnetized case does not provide any explicit graph which 
indicates that waves may not be found in the region 1.5 < z < 10. 

It would be interesting to extend this work by assuming hot plasma in 
this unusual medium. The Kerr spacetime can also be used to analyze the 
wave properties in this unusual medium. 

Appendix A 

In Veselago medium (e < 0,/i < 0), Maxwell equations, the GRMHD equa- 
tions for the general line element and the Schwarzschild planar analogue are 
given in this appendix. Maxwell equations in this unusual medium are 

V.B = 0, (Al) 

(A2) 

(A3) 

0. (A4) 
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dB 
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x E + 


~dt 


= 0, 


V 


■E = 


Je 

• 

e 


dE 


V 


x B = 


-MJ 


+ m 



The GRMHD equations in this medium will be 

f + -(B.V)H«B = --Vx(aVxB), (A5) 

dr a a 

V.B = 0, (A6) 

D Po , 2^r D V , Po f 9,t , „ ,_„ m l ft , A ^ 



Wv.^ + ^ + v.^v-^o, 

2 B 2 \ 4 1 1 D V j 

PoP7 + ^rj Hj + PoP7 ^ " ^ B i B ij -^T 

W^ - (f-^, - i-B^-) VgV* = -p o7 2 p{a 

-~iV - (iVfa)^'} -P|< + j^(V x B),V.(V x B) 

a An 

(« B )f* + ^(«^)b-^'-^:(Bx{Vx[v 



87tq; 2 /|1 ' 47ra v V|J 47TQ;' 

x (aV x B) - (B.V)/3] + (V x B).V/3}) 4 , (A8) 

n-Owy 2 ) - / + ©(PPo7 2 -P) + -(ppo7 2 ^^' 
ZJr or 2a 

+mii)£tln + 2ppo7 2 (V.a) + pp o7 2 (V.V) - V''* 

x(pp o7 V^ +P7«) - ^( V x B )-( V x ^) " ^ 

x(V x B).(B x — ) — (V x B).V/3 - — 0(V x B) 

v ' v dr 47ra v 7 4vr ' 

+— (V x O.B) = 0. (A9) 

Ana ' v ; 

Since /3, 8 and £tlij vanish for the Schwarzschild planar analogue, the perfect 
GRMHD equations reduce to 

_ = _Vx(aVxB), (A10) 

V.B = 0, (All) 

do <9V 

-g + (aV.V)po + Po7 2 V.— + po7 2 V.(«V.V)V 

+p V.(aV) = 0, (A12) 

{(PoP7 2 + 5-)*« + portWj - l-BtBM— + VS7)V* 



-^ - h BiBi)yi * yk + mVi{ ^ w + (V - VM 

= -po^ca -p,i+—(V x B)iV.(V x B) - -^(aB) 2 ;! 

47T o7rar 

+ ^ aBl),J BJ - 4^[B x {V x (V x (aV x B))}]„ (A13) 

(~ + V.V)(/i P0 7 2 ) -\% + 2m7 2 (V.a) + /i Po 7 2 (V.V) 

1 1 FfVK 1 1 r9V 

-_(V x B).(V x -^-) (V x B).(B x -^-) 

4?r V ; V a <9i ; 4vr v ; v a dt ' 

+ — (V x qB) = 0. (A14) 

Area 
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